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ABSTRACT 

We present a clustering analysis of Luminous Red Galaxies (LRGs) using nearly 9 000 ob- 
jects from the final, three year catalogue of the 2dF-SDSS LRG And QSO (2SLAQ) Sur- 
vey. We measure the redshift-space two-point correlation function, £(s) and find that, at the 
mean LRG redshift of z = 0.55, £(s) shows the characteristic downturn at small scales 
(<1 h~ 1 Mpc) expected from line-of-sight velocity dispersion. We fit a double power-law 
to £(s) and measure an amplitude and slope of s = 17.3±|;g 7i _1 Mpc, 7 = 1.03 ± 0.07 
at small scales (s < 4.5 fi^Mpc) and s Q = 9.40 ± 0.19 h^Mpc, 7 = 2.02 ± 0.07 at 
large scales (s > 4.5 /i _1 Mpc). In the semi-projected correlation function, w p (a), we find a 
simple power law with 7 = 1.83 ± 0.05 and ro = 7.30 ± 0.34 /i _1 Mpc fits the data in the 
range 0.4 < a < 50 h~ 1 Mpc, although there is evidence of a steeper power-law at smaller 
scales. A single power-law also fits the deprojected correlation function £(r), with a corre- 
lation length of ro = 7.45 ± 0.35 /i _1 Mpc and a power-law slope of 7 = 1.72 ± 0.06 in 
the 0.4 < r < 50 /i _1 Mpc range. But it is in the LRG angular correlation function that the 
strongest evidence for non-power-law features is found where a slope of 7 = —2. 17 ± 0.07 is 
seen at 1 < r < 10 /i _1 Mpc with a flatter 7 = — 1.67±0. 07 slope apparent at r < 1 ft, _1 Mpc 
scales. 

We use the simple power-law fit to the galaxy £(r), under the assumption of linear bias, to 
model the redshift space distortions in the 2-D redshift-space correlation function, £(er, tt). We 
fit for the LRG velocity dispersion, w z , the density parameter, D, m and (3(z), where (3(z) = 
r^ n 6 /6 and b is the linear bias parameter. We find values of w z — 330kms _1 , f2 m — O.lO^oio 
and (3 — 0.40 ± 0.05. The low values for w z and (3 reflect the high bias of the LRG sample. 
These high redshift results, which incorporate the Alcock-Paczynski effect and the effects of 
dynamical infall, start to break the degeneracy between fi m and (3 found in low -redshift galaxy 
surveys such as 2dFGRS. This degeneracy is further broken by introducing an additional 
external constraint, which is the value (3(z = 0.1) = 0.45 from2dFGRS, and then considering 
the evolution of clustering from z ~ to Zlrg ~ 0.55. With these combined methods we 
find n m (z = 0) = 0.30 ± 0.15 and (3(z = 0.55) = 0.45 ± 0.05. Assuming these values, we 
find a value for b(z = 0.55) = 1.66 ± 0.35. We show that this is consistent with a simple 
"high-peaks" bias prescription which assumes that LRGs have a constant co-moving density 
and their clustering evolves purely under gravity. 
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1 INTRODUCTION 

Recent measurements of the galaxy correlation function, £, have 
produced a series of impressive r esults. Whether it be th e detection 
of baryonic acoustic oscillations (Eisenstein et al. 2005! ), clustering 
prope rties of different spectral types of galaxy ( Madgwick et al 



2003), or the evolution of AGN black hole mass dCroom et al 



2005), the two-point correlation function continues to be a key 
statistic when studying galaxy clustering an d evolution. There 
have also been a ser i es of recent studies (e.g.JZehavi et alj|2005l : 
iLe Fevre et ail 120051 : ICoil et al.1 12004 IPhleps et al.l 120061) inves- 
tigating the clustering properties and evolut ion with redshift of 
galaxies from 0.3 < z < 1.5. Amongst these, Zehavi et alJ d2005l) 
use the Sloan Digital Sky Survey (SPSS; lYork et alj|200(j) to ex- 
amine the clustering properties of Luminous Red Galaxies (LRGs) 
at a redshift of z~0.35. They find that correlation length depends 
on LRG luminosity and that there is a deviation from a power-law 
in the real-space correlation function, with a dip at ~ 2 Mpc scales 
as well as an upturn on smaller scales. 

Although the form of the 2-point correlation function is in it- 
self a worthwhile cosmological datum, more information can be 
gained by studying the dynamical dist ortions at bot h small and 
large scales in the clustering pattern dKaiserl 1 19871) . Measured 
galaxy redshifts consist of a component from the Hubble expansion 
plus the motion induced by the galaxy's local potential. This leads 
to one type of distortion in redshift-space from the real-space clus- 
tering pattern. There are two basic forms of dynamical distortion 
(a) small scale virialised velocities causing elongations in redshift 
direction - 'Fingers of God', but at larger scales there will also be 
flattening of the clustering in the redshift direction due to dynami- 
cal infall. Another type of geometric distortion can be introduced if 
we assume the wrong cosmology to c onvert redshifts to comoving 
distances dAlcock & Paczvnskilll979l) . Under the assumption that 
galaxy clustering is isotropic in real-space, a test can be performed 
in redshift-space by determining which cosmological parameters 
return an isotropic clustering pattern. 

In the linear regime, dynamical effects are broadly deter- 
mined by the parameter /3, where (3 — f2 m 6 /b, fi m is the mat- 
ter density parameter and b is the linear bias parameter. If we 
assume, as is common, a zero spatial curvature model, then the 
main parameter determining geometric distortion is O m . We can 
therefore use these redshift-space distortions to our advantage and 
derive from them estimates of fi m and f3, (e.g ., Kaiser Il987t 



Lovedav et al.lll996l;lMatsubara & Sutoll996l:lMatsubara & Szalayl 
200ll;lBallinger et al.l l l996l ; |Peacock et al.ll200ll ; lHovle et al^Oolt 
da Angela et alj |2005). Unfortunately, there is often a degeneracy 
between these parameters, but this can be broken by the inclu- 
sion of other information. This additional information is introduced 
via constraints obtained fro m linear evolution the ory of cosmo- 
logical density perturbations dda Angela et al.l2005l . and references 
therein). 

In this work, we extend the redshift coverage of the SDSS 
LRG survey by using the data from the recently completed 2dF- 
SDSS LRG And QSO (2SLAQ) Survey (Cannon et al. (2006); 
Croom et al. (2007), in prep.). Luminous Red Galaxies are ideal 
candidates for galaxy redshift surveys since they are intrinsically 
bright and so can be seen to cosmological distances. Selection cri- 
teria are used which gave a relatively clean and complete selection 
of LRGs and since they are the most massive galaxies, they are 
believed to reside in over-dense peaks of the underlying matter dis- 
tribution and are thus excellent tracers of large scale- structure. 

Observations of the 2SLAQ Survey are now complete, with 



a number of new resu lts being reported (e.g. IWake et al I d2006h . 
iRoseboom e t al. (2006), Sadler et al. (2007), in prep.). In this paper 
we shall concentrate on the clustering of the 2S LAQ LRG sam- 
ple, extending the work o f the SDSS LRG Survey feisenstein et al.l 
1200 it IZehavi et alj|2005h to higher redshift. We calculate the 2- 
point galaxy correlation function in both redshift-space and real- 
space for LRGs over the redshift range 0.4 < z < 0.8. Then using 
information gained from geometric distortions in the redshift-space 



can be found (e.g. 


Alcock & Paczvnskilll979l:lBallinger et aljl996l 


IHovle et al]2002t 


da Angela et alj|2005h. 



In Section 2 we therefore introduce the 2SLAQ sample and 
the techniques used in our analysis. In Section 3 the 2SLAQ LRG 
correlation function measurements are presented and comparisons 
to other surveys are made. In Section 4 we model the redshift-space 
distortions and compare these models to our data, finding values of 
f2 m and p. Our conclusions are presented in Section 5. 



2 DATA AND TECHNIQUES 

2.1 The 2dF-SDSS LRG And QSO Survey 

A full description of the 2SLAQ Survey can be found in 
ICannon et all 0006). At its heart, the 2SLAQ Survey relies on the 
SDSS photometric survey to supply LRG targets for spectroscopic 
follow-up using the 2 Degree Field (2dF) instrument on the Anglo- 
Australian Telescope (AAT). 

The selection of distant (z > 0.4) LRGs is done on the ba- 
sis of SDSS gri photometric data, using the (g — r) versus (r — i) 
colours and the SDSS "de Vaucouleurs" i-band magnitude. The cri- 
teria are simi lar to those used for th e faint "Cut II" sample in the 
SDSS LRGs feisenstein et"aT1l200ll) and are des cribed in detail by 
ICannon et alJ d2006l) . (See iFukugita et al.l dl996T) for a description 
of the SDSS filters.) 

The survey covers two narrow stripes along the celestial equa- 
tor (\8\ < 1.5°). The Northern Stripe runs from 8^4 to 15^3 in Right 
Ascension and is broken into 5 sub-stripes to utilise the best pho- 
tometric data. The Southern Stripe runs from 20^6 to 4?0. Figure^ 
shows the layout of the target stripes and the 2dF fields observed. 
The total area of the survey, including the overlap regions, was ap- 
proximately 180 degrees 2 . Again, com plete details of the Survey 
fields are given bv lCannon et alJ |2006). 

It is important to be aware of the tiling strategy of the 2SLAQ 
survey when estimating the clustering of the LRGs. A simpler tiling 
scheme was used for 2SLAQ than for the preceding 2dFGRS/2QZ 
survey. For instance, for 2SLAQ, the 2dF tiles were offset by 1.2 
deg in the RA direction as opposed to a variable spacing strategy 
employed by the 2dFGRS and 2QZ. Again, contrary to the 2dF- 
GRS/2QZ, the galaxies in 2SLAQ were given higher fibre assign- 
ment priority, with the LRGs always having priority over the QSOs. 
This makes sure the LRG selection was not biased by the QSOs. 
The details of the survey mask and selection function will be de- 
scribed in detail in Section |23l 

The total 2SLAQ LRG dataset consists of a total of 18 487 
spectra for 14 978 discrete objects; 13 784 of these (92%) have re- 
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Figure 1. The location of the 2SLAQ Input Catalogue (dotted rectangles) and observered fields (circles). Solid circles indicate fully observed fields with high, 
> 85% overall completeness, while hollow circles have less than 85% overall completeness or fields with non-standard selection criteria. 



liable, "Qop" > 3 redshiftsQ From these "Qop"> 3 objects, 663 
are identified as being stars, leaving a total of 13 121 galaxies. 

We cut this sample down further by using only those con- 
firmed LRGs which were p art of the top priority "Sample 8" se- 
lection as described fully in ICannon et al.1 ( 2006). These galaxies 
comprise the most rigorously defined 2SLAQ LRG sample where 
completeness is highest due to their top priority for spectroscopic 
observation. The exact Sample 8 selection lines in the gri plane 
are shown in Fig. 1 of Cannon et al (2006). The magnitude limits 
is idev < 19.8 (de-reddened). However, the sample we use does 
include observations taken in the 2003A semester, where a brighter 
idev < 19.5 magnitude limit was used, as long as the observed 
LRG would have made the "Sample 8" selection. We d o not include 
observ ations taken from fields aOl, a02 and sOl (see lCannon et al.1 
(2006)) as they have low completeness and should not be used in 
statistical analyses. Once the final selection criteria had been de- 
cided, there were 25 795 "Sample 8" LRG targets at a sky den- 
sity of about 70 per square degree. Approximately 40% (10 072) of 
these objects were observed, with 9 307 obtaining "Qop"> 3. After 
imposing the cuts above, this leaves a total of 8 656 LRGs, 5 995 
in the Northern Galactic Stripe and 2 661 in the Southern Galactic 
Stripe. For all further analysis, this is the sample utilised which we 
call the "Gold Sample" and has a ZQoid = 0.55. 



Sample Description Number in sample North South 



"Qop" represents a redshift quality flag assigned by visual inspection of 
the galaxy spectrum and the redshift cross-correlation function. A value of 
3 or greater represents a 95-99% confidence that the redshift obtained from 
the spectrum is valid. 



Unique Objects 
"Qop" > 3 
M Stars 
LRGs 

LRG Sample 8 
excl. aOl, a02, sOl 



14 978 
13 784 
663 
13 121 
8 756 
8 656 



10 369 
9 726 

9 280 
6 076 
5 995 



4 609 
4 058 

3 841 
2 680 
2 661 



Table 1. The 2SLAQ LRG Survey; Numbers of different Samples. Over 
18 000 spectra were obtained, resulting in 13 121 spectroscopically con- 
firmed Luminous Red Galaxies. We use the LRGs with the "Sample 8" In- 
put Priority settings for our analysis but do not include the data taken in the 
aOl, a02 and sOl fields which have low redshift completeness and should 
be excluded from statistical analysis ( Cannon et al. 2003)- Thus we are left 
with 8 656 in our "Gold Sample". 



2.2 The Two-Point Correlation Function 

Here we give a brief description of the 2-point correlation func- 
tion (2PCF); f or a more formal treatment the reader is referred to 
IPeeblesI dl98d) which presents the basis for the rest of the section. 
To denote the redshift-space (or z-space) correlation function, we 
will use the notation £(s) and to denote the rea/-space correlation 
function, £(r) will be used, where s is the redshift-space separation 
of two galaxies and r is the real-space separation. 

The 2-point correlation function, £(x), is defined by the joint 
probability that two galaxies are found in the two volume elements 
dVi and dVi placed at separation x, 



dPn =n 2 [l + £(x)]dVidV 2 . 



(1) 



To calculate £(x), N points are given inside a window W of ob- 
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servation, which is a three-dimensional body of volume V^W). 
An estimation of t;(x) is based on an average of the counts of 
neighbours of galaxies at a given scale, or more precisely, within 
a narrow interval of sca les. An extensively used estimator is that of 
Davis & Peebles! ( fl983h and is usually called the standard estima- 
tor, 



(,std{s) = 



N rd DD{s) 
N DR(s) 



(2) 



where DD(s) is the number of pairs in a given catalogue (within 
the window W) and DR(s) is the number of pairs between the data 
and the random sample with separation in the same interval. N r d 
is the total number of random points and N is the total number of 
data points. A value of £ = 1 implies there are twice as many pairs 
of galaxies than expected for a random distribution and the scale at 
which this is the case is called the correlation length. 

2.3 Constructing a random catalogue and survey 
Completeness 

The two point correlation function, £, is measured by compar- 
ing the actual galaxy distribution to a catal ogue of randomly dis- 
trib uted galaxies. Follow ing the method of lHawkins et al.l ( 2003) 
and lRatcliffe et al.l d 1998h . these randomly distributed galaxies are 
subject to the same redshift, magnitude and mask constraints as the 
real data and we modulate the surface density of points in the ran- 
dom catalogue to follow the completeness variations. We now look 
at the various f actors this involves. 

Following ICroom et alj d2004) . we discuss issues regarding 
the 2SLAQ Survey completeness. As with the rest of the paper, we 
are only concentrating on the properties of the luminous red galax- 
ies. One might think the parallel 2SLAQ QSO survey would have 
a bearing on subsequent discussion but due to the higher priority 
given to the fibres assigned to observe the LRGs, the QSO Survey 
has no impact on LRG clustering considerations, as already noted. 
For more description of the clus tering of the QSOs the reader is 
referred to lda Angela et al. (2006). 

Three main, separate types of completeness are going to be 
considered; i) Coverage completeness, f c , which we define as the 
fraction of the Input 2SLAQ catal ogue sources that ha ve spectro- 
scopic observations. Identically to ICroom et alj d2004l) . we calcu- 
late f c , as being the ratio of observed to total sources in each of 
the sectors defined by overlapping 2SLAQ fields, which are pix- 
elized on 1 (one) arcminute scales; ii) Spectroscopic completeness, 
f s which can be said to be the fraction of observed objects which 
have a certain spectroscopic quality; iii) Incompleteness due to fi- 
bre collisions which is dealt with separately from coverage com- 
pleteness 

For coverage completeness and spectroscopic completeness 
we assume that both are functions of angular position only, i.e. 
f c {9) and f s {8) respectively. The spectroscopic (i.e. redshift) com- 
pleteness does depend on magnitude but this is not relevant for any 
of the purposes of this paper. 



2.3.1 Angular + Spectroscopic Completeness and Fibre 
Collisions 

There are various technical details associated with the 2dF instru- 
ment. Variations in target density, the small number of broken or 
otherwise unuseable fibres and constraints owing to the minimum 
fibre placing (see below) could introduce false signal into the clus- 
tering pattern. For our analysis, the 2SLAQ survey consists of 80 
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Figure 2. The w(6) for the 2SLAQ redshift catalogue (light blue) dotted, 
open circles compared to the parent catalogue solid (red) line. The errors 
quoted are "Field-to-field" errors with the sub-areas used given by Table|2] 
The filled blue squares, with dashed error bars, show the w(8) from the 
redshift catalogue after the correction for fibre collisions has been applied. 
The values for the uncorrected (corrected) w (0) from the redshift catalogue 
have been moved by A log = — (+)0.05 in the abscissa for clarity. Note 
also that the solid line represents the filled squares given in FigurefJ] 



field pointings. Many of these pointings overlap, alleviating some 
of these technical issues. 

The design of the 2dF instrument mean s that fibres cannot be 
placed closer than approximately 30 arcsec dLewis et al.ll2002h so 
both members of a close pair of galaxies cannot be targeted in a 
single fibre configuration. The simple, fixed-spacing tiling strategy 
of the 2SLAQ Survey means that not all such close pairs are lost. 
Neighbouring tiles have significant areas of overlap and much of 
the survey sky area is targeted more than once. This allows us to 
target both galaxies in some close pairs. Nevertheless, the survey 
misses a noticeable fraction of close pairs. It is important to assess 
the impact of this omission on the measurement of galaxy cluster- 
ing and to investigate schemes that can compensate for the loss of 
close pairs. 

To quantify the effect of these so-cal led 'fibre collisions' 
we have followed previous 2dF studies (e.g. lHawkins et al]|2003l ; 
ICroom et ai]|2004l) and calculated the angular correlation function 
for galaxies in the 2SLAQ parent catalogue, w p (6), and for galax- 
ies with redshifts used in our £ analysis, w z {9). We used the same 
mask to determine the angular selection for each sample. 

As shown in Figure [2] on scales 9 > 3', the angular correla- 
tions of the Parent and Redshift catalogue are very nearly consis- 
tent. At scales 8 < 2', we begin to lose close pairs. To corre ct for 
this effect, we use a similar method to Hawki nTet alj d2003l) and 
iLi et alj J2006l) . The quantity w COT (6) = (1 + w p )/(l + w z ) is 
used to weight our 3-D DD pairs. For each DD pair, the angular 
separation on the sky is calculated and the galaxy-galaxy pair is 
weighted by the w COI (0) ratio given by the relevant angular sepa- 
ration. The result of weighting by this factor, is shown by the filled 
(dark blue) squares in Fig. [2] 

The last stage in determining the angular "mask" is to evalu- 
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Figure 3. The redshift distribution for the 2SLAQ LRG "Gold" Sample we 
use. The solid red histogram is for the "Gold" Sample. The dashed blue line 
is from the normalised random catalogue. 



ate the spectroscopic completeness of the survey, f s (8) which for 
our purposes, we again assume depends on sky position only. This 
function essentially describes the success rate in obtaining a spec- 
trum and reliable redshift for a given fibred object. Here the ad- 
vantage of LRGs becomes apparent. With their well-defined early- 
type spectra and often very strong Ca H+K break around 4000A, 
a high success rate was achieved when calculating a redshift for 
the 2SLAQ LRG objects. Also, it became apparent that our 4 hour 
per field exposure time was on occasion generous and relatively 
high S/N spectra were recorded. The spectroscopic completeness 
has been estimated at 94.5 per cent for the primary "Sample 8" and 
the redshift completeness at 96.7 per cent, giving an overall com- 
pleteness of 91.4 per cent (Cannon et al. 2006, Section 5.5, Figure 
5). 



2.3.2 Radial Selection Function and Estimates of the LRG N(z) 

The observed distribution of galaxy redshifts is given in Figure [3] 
Plotted are the N(z) distributions, binned into redshift slices of 
Az=0.02, for the "Gold Sample". Also shown is a polynomial fit 
(7th order) to the N(z) distribution, which is used to generate the 
random distributions. Checking the N(z) fits using higher order 
polynomials or convolved double Gaussians does not give tighter 
reproduction of the observed LRG redshift distribution. 

Combining the radial selection function and the completeness 
map, we generate a random catalogue of points which we now use 
to calculate the LRG correlation function. 



2.4 Calculating the 2-point Correlation Function 

As the LRG correlation function, £(s), probes high redshifts and 
large scales, the measured values are highly dependent on the as- 
sumed cosmology. In determining the comoving separation of pairs 
of LRGs we choose to calculate £(s) for two representative cos- 
mological models. The first uses the cosmolo gical parameters de- 
rived from WMAP, 2dFGRS and other data dSpergel et ai]|2003l 



l2006l : |Percivaletal]|2002l ; ICole et al.l [20051 : ISanchez etai]|2006l) 
with (fl m ,fiA) — (0.3, 0.7), which we will call the A cosmology. 
The second model assumed is an Einstein-de Sitter cosmology with 
(fi m ,S7A) = (1.0, 0.0) which we denote as the EdS cosmology. We 
will quote distances in terms of h^ 1 Mpc, where h is the dimen- 
sionless Hubble constant such that Ho = 100/t km s~ Mpc - . 

We have used the m inimum variance estimator suggested 
by Landv & Sza lavl (11993b to calculate £(s). Using notation from 
Marti nez & Saarj ( 120021) . this estimator is 
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where the angle brackets denote the suitably normalised LRG- 
LRG, LRG-random and random-random pairs counted at separa- 
tion s. We use bin widths of <51og(s/ /i -1 Mpc) = 0.1. The den- 
sity of random points used was 20 times the densi ty of LRGs. The 
Hamilton estimator is also utilised dHamiltonll 19931) where 

DD(s)-RR(8) 

5Ham(Sj — nnr „ 1 (5) 



DR(s) 2 

and no normalisation is required. Since we find the differences of 
the Hamilton estimator compared to the Landy-Szalay method are 
negligible, the Landy-Szalay method is quoted in all £(s) figures 
unless explicitly stated otherwise. 

Three methods are employed to estimate the likely errors on 
our measurements. The first is a calculation of the error on £(s) 
using the Poisson estimate of 



CPo 



!+€(«) 



(6) 



The second error estimate method is what we shall call the 
field-to-field errors, calculated by 



N - 



DRjjs) 
DR(s) 



(7) 



where N is the total number of subsamples i.e. "the fields" and 
£i(s) is from one field. { (s) is the value for £ from the entire sample 
and is not the mean of the subsamples. For our studies the natural 
unit of the "Field-to-field" (FtF) subsample is given by the area 
geometry covered by the survey. Thus we take N — 9, and split the 
NGP area into five regions, a,b,c,d,e and the SGP in to four regions, 
named s06, s25, s48, s67. Details of the FtF subsamples are given 
in Table|2 

The third method is usually referred to as the jackknife 
estimate, and has been used in other corre lation studies (e.g. 
IScranton et aT]|2002l : liehavi et al.ll200i l2005h . Here we estimate 



fjacklS) 
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(8) 



where i' is used to signify the fact that each time we calculate a 
value of £(s), all subsamples are used bar one. For the jackknife 
errors, we divide the survey into 32 approximately equal sized ar- 
eas, leaving out ~4.5 square degrees from the entire survey area at 
one time. Thus a jackknife subsample will contain ~8,350 LRGs. 
We can then work out the covariance matrix in the traditional way, 
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(9) 
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Figure 4. The ratio of Poisson to jackknife errors (solid black line and 
squares), Poisson to 'field-to-field' errors, (dashed blue line and triangles) 
and 'field-to-field' to jackknife errors (dotted red line and open squares). As 
can be seen, all error estimators are comparable on scales < 10 h~ 1 Mpc, 
while at larger scales than this the jackknife and 'field-to-field' errors are 
considerably larger than the simple Poisson estimates. The magnitude of 
the 'field-to-field' and jackknife errors are very similar from the smallest 
scales considered here up to Rj 40 /i _1 Mpc. 



where £ is the mean value of £ measu red from all the jac kknife 
subsamples and TV = 32 in our case (c.f. lZehavi etail J2002h ). The 
variances are obtained from the leading diagonal elements of the 
covariance matrix, 



CFi = COV 



(10) 



When examing the covariance matrix, we find the measurements 
to be slightly noisy as well as an indication of anti-correlation 
(contrary to theoretical expectations). However, we note that in 
the other recent clustering stu dies, noisy covariances and anti- 
corre l ations were also noted (e.g. lScranton et al"]2 002; Zeha vi et al.l 
l200ll2005h . 

The ratio of Poisson to jackknife errors, Poisson to 'field-to- 
field' errors, and the 'field-to-field' to jackknife errors are given in 
Figure [4] As can be seen, all error estimators are comparable on 
scales < 10 /i _1 Mpc, while at larger scales than this the jackknife 
and 'field-to-field' errors are considerably larger than the simple 
Poisson estimates. The magnitude of the 'field-to-field' and jack- 
knife errors are very similar from the smallest scales considered 
here, up to w 40 /i -1 Mpc. This b ehaviour has been no ted in 
other correlation function work, e.g. Ida Angela et al. | fe005h . We 
also note that field-to-field and jackknife errors are more compa- 
rable in size, regardless of scale. Hence, the errors that are quoted 
on all correlation functions from here on are the square roots of 
the variances from the jackknife method, except for the case of the 
angular correlation function, w(8), where we quote the "Field-to- 
field" error. 



Area Name 


RA(J2()00) range/ 


LRGs 


Randoms 


O J 1 flT T>/~i 

Hra/ yLttikj 


a 


123.0- 144.0 


617 


10 745 


17.41 


b 


150.0- 168.0 


1 837 


35 449 


19.30 


c 


185.0- 193.0 


572 


14 484 


25.32 


d 


197.0-214.0 


1 723 


34 373 


19.95 


e 


218.0- 230.0 


1 246 


24 849 


19.94 


s()6 


309.2 - 330.0 


745 


12 457 


16.72 


s25 


330.0 - 360.0 


876 


18 499 


21.12 


s48 


0.0 - 30.0 


658 


13 516 


20.54 


s67 


30.0 - 59.7 


382 


8 749 


22.90 


Entire Survey 




8 656 


173 120 


20.00 



Table 2. The 2SLAQ LRG Survey; Names and Right Ascension ranges for 
the N = 9 sections used when calculating the field-to-field errors. 



2.5 Measuring £(<r, tt) 

Having described how we calculate galaxy-galaxy separations in 
redshift-space in order to measure f (s), we can now study the clus- 
tering perpendicular, a, and parallel, tt, to the line of sight. We work 
out the co-moving distance, r c , to our object, which is equal to the 
distance parallel to the line of sight i.e. a tt value. Thus, already 
knowing the redshift-space separation, s, we can use 



2 2,2 
s — a + TT 



(11) 



to find a. At this point it should be noted that a is sometimes 
designated by r p , where r p = a. For this paper we shall con- 
tinue to use a for the perpendicular separation. Closely following 
iHovle et alj d2002h . £(<r, tt) can be estimated in a similar way to 
£(s). A catalogue of points, that have the same radial selection 
function and angular mask as the data but are unciustered, is used 
to estimate the effective volume of each bin. As stated above, the 
unciustered, random catalogue also contains 20 times more points 
than the data. The DD(a, tt), DR(a, tt) and the RR(a, tt), where 
again D stands for data LRG and R stands for random, counts in 
each a and tt bins are found and the Landy-Szalay estimator 



£ls{<t, tt) 



{DD(tj, tt)} - (2DR(o, tt)) + {RR{a, tt)) 
(RR(a,Tr)) 



(12) 



is used to find f(cr,7r), with bins of 5\og(a/ /i _1 Mpc) = 
<51og(7r/ /i~ 1 Mpc) = 0.2. Again, we compute three types of er- 
rors to use as a guide; Poisson, "Field-to-field" and Jackknife er- 
rors are calculated for £(<r, tt) as in equations [6] to [8] Again, after 
comparing the different £(cr, tt) error estimators we find that on the 
scales we are considering, the jackknife error is sufficient for our 
purposes. 



2.6 The Projected Correlation Function, w p (a) 

Although we are now in a position to calculate the redshift-space 
correlation function, the real-space correlation function, £(r), 
which measures the physical clustering of galaxies and is indepen- 
dent of redshift-space distortions, remains unknown. However, due 
to the fact that redshift distortion effects only appear in the radial 
component, by integrating along the n direction, we can calculate 
the projected correlation function, 



w p (a) 



£(cr, 7r)d7r 



(13) 
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In practice we set the upper limit on the integral to be 7r max = 
70 Mpc as at this large-scale, the effect of clustering is negli- 
gible, while linear theory should also apply. The effect of z-space 
distortions due to small-scale peculiar velocities or redshift errors 
is also minimal on this scale. Changing the value of 7r max from 
25 Mpc to 100 h" 1 Mpc makes negligible difference in the 
result. 

Due to Wp(a) now describing the real-space clustering, 
the integral in Equatio n [T3l can be re-written in terms of f (r), 
jDavis & Peeblesll 19831) 



w p (a) 



J a 



r£(r) 



zdr. 



(14) 



If we then assume that £ (r) is a power-law of the form, £(r) = 
( r / r o) 1 , equation[l4]can be integrated analytically such that 



w p (a) 



r(|)r( 



-y — 1 * 



(15) 



where A(pj) represents the quantity inside the square brackets and 
T(x) is the Gamma function calculated at x. We now have a method 
for calculating the real-space correlation length and power-law 
slope, denoted ro and 7 respectively. 



2.7 The Real-space Correlation Function, f (r) 

Using the projected correlation function, w p (a), it is now possi- 
ble to find the ro and 7 for the real-space correlation function. 
However, if one does not assume a power-law f(r), it is still 
possible to estimate £(r) by directly inverting w p (a). Following 
ISaundersetafldl992l) we can write 



f(r) = -r 



(dw(a) /da) 



da. 



(16) 



Assuming a step function for w p (a) — Wi in bins centred on at 
and interpolating between values, 



■7T ' ■* 



Wj+i 



In 



cr + 



(17) 



for r = ai. We shall be utilising this interpolation method to check 
whether a power-law description is valid for our 2SLAQ Survey 
data and, if so, what values the parameters ro and 7 take. 



3 RESULTS 

3.1 The LRG Angular Correlation Function, w(9) 

We first analyse the form of the angular correlation function, w(9). 
The full input catalogue contains approximately 75 000 LRGs 
mainly from areas in the two equatorial stripes; about 40% of this 
area was observed spectroscopically. As stated in Section 2, ap- 
proximately a third of the objects in the total input catalogue pass 
the Sample 8 selection criteria. As well as providing estimates of fi- 
bre collision and other angular incompletenesses, the angular func- 
tion is of interest in itself, particularly given the narrow redshift 
range from which the sample is derived. We use 25 795 "Sample 
8" LRG targets to estimate the w(9). We first note that the func- 
tion gives clear indication of a change of slope at 9 = 2 arcmin 
or ~ 1 ft - 1 Mpc in the A cosmology. Considering the form of 
w{9) = A9 1 ~~<, at 9 < 2 arcmin the slope is 7 = -2.17 ± 0.07 



and at larger scales the slope is 7 = —1. 67 ± 0.07. Using Lim- 
ber's formula from Philli pps et al] Jl978h and assuming a dou- 
ble power-law form where the slope changed between -2.17 and 
-1.67 at ~ 1.5 h~ Mpc (comoving), we found in the A case, a 
value of ro = 4.85 ± 0.3 /i -1 Mpc at small scales and ro = 
6.89 ± 0.6 h~ Mpc at large scales (see Fig. [5]l. We shall check 
models of this form against the deprojected correlation function 
£(r) (see Figure [9] below). We find that the form of this double 
power-law gives reasonable fits to the data in the LRG redshift sur- 
vey, although the large scale slope derived from the input catalogue 
w(9) appears slightly flatter than in the semi-projected and 3-D 
correlation functions (see below). The reason for this is not clear, 
although it could be that w[ff) is more sensitive to any artificial 
gradient in the LRG data. Thus, we checked for an angular sys- 
tematic in the data by calculating the angular correlation between 
spectroscopic LRGs that are not at the same redshift. We find this is 
consistent with zero and so such systematics do not explain the flat- 
ter slope for w(9) at large-scales. The most likely explanation is the 
different fitting ranges for w(9) and the semi-projected correlation 
function. This test also suggests that the upturn at 9 < 2 arcmins 
is a real feature. It will be seen that w(9) gives the strongest ev- 
idence of all the correlation function statistics for non-power-law 
behaviour in £ (r). A similar feature is seen by Zehavi et al in the 
SDSS MAIN galaxy sample and to a lesser extent in the SDSS LRG 
survey. Reports of such fea tures in galaxy correlation functions go 
back to lShanks et al. I dl983h . We simply report the existence of this 
feature in the LRG data and leave further interpretation for a fu- 
ture paper. Possible interpretations could include models of halo 
occupation distributions (HOD) in the standard model case or the 
possibility that it might represent a real feature in the mass dis- 
tribution in the ca se of other models. We also show results from 
IWhite et al.l 12001 open, black circles, Figure[5]( who report on the 
angular correlation function as a route to estimating merger rates of 
massive red galaxies. As can be seen, these measurements from the 
NOAO Deep Wide-Field Survey (NDWFS: |jannuzi & Devlll999l) 
agree very well with the 2SLAQ LRG results, though as we shall 
discuss later, care always has to be taken when comparing measure- 
ments from galaxy surveys with different selections. 

3.2 The LRG Redshift-Space Correlation Function, £(s) 

Using the above corrections including that for fibre collisions, 
the 2SLAQ LRG redshift-space 2PCF, £(s), is shown in Fig- 
ure [6] There is clear evidence for downturns at small scales 
<2.5 /i _1 Mpc and large scales >10 /i _1 Mpc that are not de- 
scribed well by a single power-law. This turn-over is consistent with 
the redshift-space distortion effects one would expect in a f (s) cor- 
relation function - namely the "Finger of God" effect at small scales 
due to intrinsic velocity dispersions and large-scale flattening from 
peculiar motions due to coherent cluster in-fall. However, we note 
that real features in the real-space correlation function, f (r), may 
also be contributin g. We have als o estimated the effect of the inte- 
gral constraint (IC, Peebles 1980) at larger scales. Using our global 
(N+S) normalisation of the correlation function, we assume a total 
number of 8 656 galaxies in a total volume of 4.5 x 10 7 /i -1 Mpc 3 
and ro = 7.45 /i -1 Mpc. Integrating with a 7 = 1.8 power-law to 
20 fc _1 Mpc gives an IC = 3.5 x 10" 4 and to 100 /i _1 Mpc, an 
IC — 2.4 x 10~ 3 . Adding such contributions would make negli- 
gible contributions to any of our correlation function fits. 

We now attempt to parameterise the £(s) data. The simplest 
model traditionally fitted to correlation function estimates is a 
power law of the form 
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Figure 5. The angular correlation function, w(8) from the 2SLAQ input 
catalogue containing 25 795 LRG targets (solid, red squares). Clear evi- 
dence is seen for a change of power-law slope on ~ 2arcmin scales which 
is equivalent to 1 /i^M pc. The open (b lack) circles show the results 
from the NDWFS at z ~ 0.5 Iwhite et alj2007l) . 
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Figure 7. The redshift-space correlation function, f (s) for the 2SLAQ LRG 
Survey (filled, red, diamonds). For comparison, data from th e SDSS LRG 
Survey (black stars Zehavi et al. 2005; Eisenstein et al. 2005) and the high 
luminosity early-type 2dFGRS, INorberg et al.ll2002l open blue triangles) 
are also plotted. 
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Figure 6. The redshift-space 2-point correlation function, f(s) for the 
2SLAQ LRG Survey in a A cosmology (filled, red diamonds) and an 
Einstein-de Sitter, n m = 1, cosmology (open, cyan diamonds). The dashed 
lines shown are the double power-law best-fit models to data with the asso- 
ciated values of so and 7 given in Table|5] 



(18) 



where so is the comoving correlation length, in units of h~ 1 Mpc. 
However, with the redshift-space distortion effects being so evi- 
dent, we find that a single-power is insufficient to describe the data 
and thus switch to a double power-law model 



Table 3. Values of the redshift-space correlation length and slope for the 
2SLAQ LRG Survey from £(s). When a A cosmology was assumed, s b 
was set at 4.5 /i _1 Mpc. When a EdS cosmology was assumed, s;, was set 
at 2.5 /i _1 Mpc. 



A 




so < 4.5 /i _1 Mpc 


so > 4.5 h _1 Mpc 


s / h 

7 

x 2 ■ 


_1 Mpc 
(reduced) 


I7.3t|; 
1.03 ±0.07 
1.95 


9.40 ±0.19 
2.02 ±0.07 
1.88 


EdS 




so < 2.5 /i _1 Mpc 


s > 2.5 h _1 Mpc 


so/ h 

1 

x 2 

^rmn 


^Mpc 
(reduced) 


20.3t|; 
0.88 ±0.11 
0.91 


7.15 ±0.13 

1 qq+O-05 

3.43 



s < Sh and 



s > s h 



(19) 



where Sb is the scale of the "break" from one power-law descrip- 
tion to the other. This f(s) model is used later in Section [4] We 
fit the double power-law continuously over the range 0.4 < s < 
60 /i _1 Mpc. We fix the break-scale at 4.5 /i _1 Mpc for the A cos- 
mology and at 2.5 /i _1 Mpc for the EdS c osmology. We per form a 
X 2 -fit, following the prescription given bv lPress et al.l (l992, Chap. 
15), to find the best-fit values for si,7i,S2, and 72. We plot the 
best fit double-power law models in Figure [6] and quote the val- 
ues of si, 71, s%, and 72, in Table|3] The errors quoted in Table 3 
are only indicative because no account has been taken of the non- 
independence of the correlation function points in deriving the £(s) 
fits. 

For com parison, in Figure[7]results from the SDSS LRG study 
are reported dZehavi et alj 20051 : iEisenstein et al. 2005) as well as 
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Table 4. Values of the projected correlation function, w p (cr), correlation 
length and slope for the 2SLAQ LRG Survey. In the A model, fits were 
performed over the range 0.4 < cr < 50.0 /i^Mpc, whereas for the EdS 
model, fits were performed over 0.25 < a < 40.0 /i _1 Mpc. The value of 
rg was found using equation [T5l 





A 


EdS 


r / h~ 1 Mpc 


7.30 ±0.34 


5.40 ±0.31 


7 


1.83 ±0.05 


1.82 ±0.06 


Xmin (reduced) 


1.17 


1.39 



selected measurements from the 2dFGRS dNorberg et ail |2002|) . 
The 2dFGRS is a blue, b.j selected sur vey of generally ~ L* galax- 
ies. However, in lNorberg et alj l l2002h . the sample is segregated by 
luminosity and spectral type, the latter governed by the r\ parameter 
(Madg wicketal.l2003h . Assuming a conversion of M° 2 -M b , ~ 
— 1.1, we calculate that the faintest 2SLAQ LRGs in our sample 
have a n M bl ~ —20 . 5. We ighting according to number, we thus 
use the lNorberg et alj |20o3) -21.00 > M bj - 5 log ft > -22.00 
and —20.50 > Mbj — 5 log h > —21.50 luminosity ranges from 
their "early-type" volume-limited sample. This is shown by the 
(blue) open triangles in Figure[7] 

The 2SLAQ LRG measurement is lower than the SDSS LRG 
result. It should not be concluded that this is evidence of evolu- 
tion because although the SDSS survey is at a lower mean red- 
shift, it was desi gned in order to targe t generally redder, more lu- 
minous LRGs dEisenstein et alj 12001). The 2SLAQ LRG colour 
selection criteria is relatively relaxed for an "LRG" survey, lead- 
ing to bluer and less luminous galaxies making it into our sam- 
ple. We note here that it is non-trival comparing clustering ampli- 
tudes and bias strengths for surveys with (sometimes very) differ- 
ent colour/magnitude/redshift selections. As such, a more detailed 
analysis of the clustering evolution for SDSS and 2SLAQ LRGs is 
presented in Wake et al. (2007, in prep.). 

The 2dFGRS M b] < —20.5, early-type sample is at least ap- 
proximately matched in terms of luminosity to the 2SLAQ LRGs. 
Once we have determined the linear bias parameter b for the z = 
0.55 2SLAQ LRGs, we shall be able to use a simple model of bias 
evolution, to compare these low redshift 2dFGRS and 2SLAQ LRG 
results. 



3.3 The Projected Correlation function, w p (a) 

Again, after applying coverage, spectroscopic and fibre collision 
corrections, the projected correlation function, w p (a), is presented 
in Figure [8] We again fit a single power-law to the 2SLAQ data 
and find that for the A cosmology, a single power-law is an ade- 
quate description, returning a reduced \ 2 = 1-17 over 0.4 < a < 
50 ft _1 Mpc. Over the wider range of 0.1 < a < 50 h~ 1 Mpc, 
the x 2 increases to 1.71. Thus the projected correlation function 
appears to deviate from a single power law at small scales in the 
way described in Section [3~T] The results for ro and 7 assuming a 
single power-law are given in Table [4] The errors are taken from 
jack-knife estimates found by dividing the survey into 32 subareas. 

This power-law deviation in the projected correlation func- 
tion is in line with recent results seen in oth er galaxy surveys, e.g. 
the SDSS M AIN sample dZehavi et alj d2004h . not plotted) and the 
SDSS LRGs dZehavi et al|2005» . A ' 'shoulder" is reported in these 
studies around ~ 1 ft _1 Mpc scales. This feature is currently be- 
lieved to be a consequence of the transition from the measuring of 
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Figure 8. The 2SLAQ LRG projected correlation function, w p (cr), with er- 
rorbars from the "Jackknife" estimates (solid, red diamonds). The dashed 
line is the power-law that gives the best fitting line from t he x 2 analy- 
sis (s ee Table [4}- The measurements from the SDSS LRGs (Zehavi et al. 
120051) are shown as a guide, with the SDSS errors being of comparable 
size to the plo tted stars. The op en (green) triangles are from COMBO- 17 
Red Sequence (Phl eps et al J2006I) . The lower panel shows the 2SLAQ LRG 
w p (a) measurements divided by this best- fitting power law with the dashed 
line covering 0.4 < a < 50 h~ Mpc. 



galaxies that reside within the same halo (the "one-halo" term) to 
the measuring of galaxies in separate haloes (the "two-halo" term). 
Dips in the projected correlation function are a major prediction of 
HOD models. Thus for the 2SLAQ LRG Survey, we set a fiducial 
model, based on our best-fitting single power-law model of w p (a) 
and find that if we divide the data out by this model, th e results (bot- 
tom pa nel, Figure[8j are potentially comparable to the lZehavi et al.1 
(2005) results (their Figure 1 1). Despite the fact that our LRG sam- 
ple is at higher redshifts and extends to lower luminosities, the form 
of the projected correlation function appears close to that seen in 
the SDSS LRG sample, although at lower amplitude. We conclude 
that the 2SLAQ LRG correlation function changes slope in similar 
fashion to the SDSS LRG semi-projected correlation function. 

Continuing with w p (a), we compare the 2SLAQ LRGs with 
the COMBO- 17 Survey. C OMBO-17 (Clas sifying Objects by 
Medium-Band Observations, IWolf et al J J200lh uses a combination 
of 17 filters to obtain photometric redshifts accurate to <r z /(l + 
z) ~ 0.01 for the brightest (i?v cg a < 20 mag) objects. This is a 
comparable sample to our own in that it covers the same redshift 
range (0.4 < z < 0.8), but care must be taken when compar- 
ing the results; although the COMBO-17 galaxies described here 
are defined as Red Sequence, on the whole they will not be LRGs 
and will have a fainter magnitude and different colour selection. 
Figure [8] gives the projected correlation function of t he 2S LAQ 
LRGs and red COMBO-17 galaxies from lPhleps et alj d2006t) (as- 
suming a flat A cosmology). The change in slope is clearly seen 
in COMBO- 17 and indeed is m odelled successfully with a HOD 
prescription dPhleps et alJl2006h . The upturn in slope in COMBO- 
17 versus 2SLAQ seems to occur on slightly different scales (~ 
1 — 2 /i _1 Mpc versus ~ 5 ft _1 Mpc) and is more dramatic than 
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Figure 9. The real-space 2-point correlation function for the 2SLAQ LRG 
Survey (filled, red, diamonds) for the A cosmology. The best-fit single 
power-law with ro = 7.45 ± 0.35 and 7 = 1.72 ± 0.06 is given by 
the dashed (red) line. The double power-law fit reported for the angular 
correlation, w(8), in Section 3.1, is shown by the dotted (blue) line. The 
solid (black) line is a theoretical predict ion for the £mass(.z = 0.55) using 
the simulations from lCoh'n et alj Il999l) . These models have (n m ,f2A) = 
(0.3, 0.7), h = 0.7 and a as = 1.0. We shall return to this in Section 4. 
The lower panel shows the 2SLAQ LRG £(r) measurements (assuming a 
A cosmology) divided by this best-fitting power law with the dashed line 
covering 0.4 < a < 50 /i~ 1 Mpc. 



Table 5. Values of the correlation length and slope for the 2SLAQ LRG 
Survey from the real-space correlation function, £(r). Model fits were per- 
formed over the range 0.4 < r < 50 h _1 Mpc for the A cosmology and 
over the range 0.25 < r < 40 h~ 1 Mpc for the EdS cosmology. 





A 


EdS 


r / h _1 Mpc 


7.45 ±0.35 


5.65 ±0.41 


7 


1.72 ±0.06 


1.67 ±0.09 


Y 2 . (reduced) 


1.73 


0.62 



for either of the LRG samples. The errors on the COMBO- 17 data 
are also much greater. Whether the differences are real, caused by 
the fainter magnitude of the COMBO-17 galaxies, or whether they 
are due to anomalies caused by the photometric redshifts, remains 
unclear. 

3.4 The Real-Space Correlation function, f (r) 

Having reported the clustering of 2SLAQ LRGs using the z-space 
correlation function, f(s) and the projected correlation function, 
w p (a) we now use the methods quoted in Section 2 to estimate the 
real-space correlation function, £(r), We show this in Figure|9] 

Again, we attempt to fit simple power-law models to our £(r) 
data in order to find values for the real-space correlation length 
and slope, ro and 7, respectively. For £(r) we attempt to take into 
account the information presented in the covariance matrix by esti- 
mating x 2 fits t° model £(r) values such that 



i,j 

where CZ l is the inverse matrix of the covariance matrix and the 
subscripts i and j are indicies of separation bins. Ho wever, as has 
been reported in prev i ous c lustering analyses (e.g. IZehavi et al.l 
d2002l) : IScranton et alj {200$) ), the calculated covariance matrix is 
rather noisy with anti-correlations between points (contary to the- 
oretical expectations). Therefore, when calculating the best-fitting 
models, we perform a simple \ 2 fit as before, without the covari- 
ances, and take only the variances into account. As before, we fit 
over the scales 0.4 < r < 50.0 h~ Mpc. For the case of the real- 
space correlation function, we again find that a single power-law 
may not fit the data well with the best-fit values (and related re- 
duced \ 2 ) given in Table|5] We find a value of 7 to be 1.72 ± 0.06 
and a correlation length of ro = 7.45±0.35 (assuming a A cosmol- 
ogy). The errors on these parameters are estimated from consider- 
ing the la deivation from the minimized \ 2 on the 1-parameter fits. 
However, care has to be taken when quoting the best fit values for 
the joint 2-parameter fits which are shown in Figure [10] Here we 
find the values of 8\ 2 which correspond to the 1, 2 and 3a levels for 
a 2-parameter fit. Also shown in Fig. [10] are the values for the de- 
viations in ro and 7, if we find the 32 best-fitting single power-law 
parameters from the jackknife samples. Jackknife appears to con- 
firm the \ 2 error analysis with the assumption of Gaussian errors 
in Fig. [TO] This is somewhat surprising since we have ignored the 
covariance between correlation function points in creating Fig. [10] 
The explanation may be that the fit at the minimum is still poor due 
to the deviant point at 2 ft _1 Mpc in Fig.|9]and this causes the error 
contours in Fig.[H)]to be larger than they would be in the absence 
of the deviant point. Including the full covariance matrix, the A^ 2 
produces error contours significantly smaller than those in Fig. 1 101 
and also the jackknife errors, even though the \ 2 at minimum re- 
mained the same. Overall we take the errors in Fig. 1 101 supported 
by the jackknife estimates as being reasonably representative of the 
real error. 

Now armed with our best-fitting single power-law model for 
£(r), and we can proceed and see if modelling the redshift-space 
distortions introduced into the clustering pattern reveals anything 
about cosmological parameters. 



4 LRG CLUSTERING AND COSMOLOGICAL 
IMPLICATIONS 

Having calculated the z-space, projected and real-space correlation 
functions for the 2SLAQ Luminous Red Galaxies, we can now turn 
our attention to using these results to see if we can determine cos- 
mological parameters. 

4.1 The £(cr, 7r) LRG Measurements 

Results for the 2-D clustering of 2SLAQ LRGs are shown in the 
£(<7, 7r) plots of Figures 1 1 1 land[T2l 

Galaxy peculiar velocities lead to distortions in the n) 
shape. The predominant effect at large scales in a is the coherent in- 
fall that causes a flattening of the £(cr, 71") contours along the parallel 
n direction and some elongation along the perpendicular a direc- 
tion. At small a, the random peculiar motions of the galaxies cause 
an elongation of the clustering signal along the 7r direction - the so- 
called "Fingers-of-God" effect. From the measurements of these 
effects, determination of the coherent infall into clusters, given by 
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Figure 10. The joint 2 parameter fits on ro and 7 for £(r). The contours 
show the 8\ 2 = (2.3,6.17,11.8) corresponding to 1, 2 and 3cr. The 
crosses show the deviations in ro and 7 that we find from the 32 best-fitting 
single power-law using the jackknife samples. 
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Figure 12. The f (ct, tt) contour plot for the 2SLAQ LRG Survey, with a 
Qui = 1.0, EdS cosmology. 

{w 2 ) respectively. We then use f3 to find the bias of LRGs at the 
survey redshift. Next, we employ information gained in studying 
the geometric distortions to perform the "Alcock-Paczynski Test" 
as one route to calculating cosmological parameters. However, we 
realise there is a degeneracy in the (/3, fi m ) plane with this ap- 
proach and thus employ further constraints from the evolution of 
LRG clustering to break this degeneracy. 




Figure 11. The £ (a, it) contour plot for the 2SLAQ LRG Survey, assuming 
a A cosmology of (f2 m , Qa) = (0.3, 0.7). The "Finger-of-God" effects, 
i.e. elongation of contours in the tt direction at small (<1 h _1 Mpc) scales, 
are clearly seen. (The spikes at small tt are a plotting artefacts). 



the parameter (3, and the pairwise velocity dispersion, (w 2 ) 1 , can 



be made. This calculation shall be performed in Section |472| Geo- 
metric distortions also occur if the cosmology assumed to convert 
the observed galaxy redshifts is not the same as the true, underlying 
cosmology of the Universe. The reason for this is because the cos- 
mology dependence of the separations along the redshift direction 
is not the sam e as for the separations me asured in the perpendic- 
ular direction dAlcock & Paczvnski| [l979V We note that modelling 
the geometric distortions and comparing to the presented data can 
yield information on cosmological parameters. 

We shall closely follow the methods of Hoyle et al. (2002) 
and da Angela (2005), hereafter H02 and dA05, respectively. In 
this section, we first discuss large-scale, linear and small-scale non- 
linear z-space distortions and how they are parameterised by and 



4.2 Redshift-space distortions, j3 and pairwise velocities 

When measuring a galaxy redshift, one is actually measuring a sum 
of velocites0 The total velocity comes from the Hubble expansion 
plus the motion induced by the galaxy's local potential, where this 
second term is coined the "peculiar velocity", i.e. 



«Tot = VH + V pec 

The peculiar velocity itself contains two terms, 

Vpec = U r and + W CI 



(21) 



(22) 



The first term, v lan d is due to the small-scale random motion of 
galaxies within clusters. The second term, vci is the component 
due to coherent infall around clusters, where the infall is caused by 
the streaming of matter from underdense to overdense regions; this 
leads to a "flattening" in the perpendicular cr-direction away from 
equi-distant contours in £(er, tt). This extension is parameterised by 
(5, which takes into account the large-scale effects of linear z-space 
distortions. Kaiser (1987) showed that, assuming a pure power-law 
model for the real-space correlation function (which is fair for the 
2SLAQ LRG data), one can estimate j3 in the linear regime using 



£(*) = £(r) 



2 1 , 



(23) 



and more generally 
2(1 



£(a )7 r) = 



1 



7 /i 2 )„ , 3-6 7 /i 2 +7(2 + 7)^ 



(3- 7 )(5-7) 



£(0,(24) 



where /i is the cosine of the angle between r and tt (the dis- 
tance along the line of sight), and 7 is slope of the power law 
jMatsubara & Sutolll996t) . 



2 This section strongly follows Hawkins et al. 12003) and Croom et al. 
120051) . 
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Even though the "Kaiser Limit" is a widely used method for 
estimating 0, drawbacks using this approach, under the assumption 
of Ga ussianity, have been known for some time dHatton & Cold 
1998). Scoccimarro (2004) has recently reported on the limitations 
of assuming a Gaussian distribution in the pairwise velocity disper- 
sion <7i2, even at very large scales. Scoccimarro's argument is that 
even at large scales, linear theory cannot be applied since one still 
has the effect of galactic motions induced on sub-halo scales i.e. 
galaxies that are separated by very large distances are still "hum- 
ming" about inside their own dark mattter haloes. Thus for the re- 
mainder of the paper, we make a note of the new formalism in Scoc- 
cimarro (2004), but continue to use the Kaiser Limit, acknowledg- 
ing its short-comings. We justify this by noting that we need better 
control on our '1st order' statistical and systematic errors before 
applying the '2nd order' Scoccimarro corrections. Future analysis 
may use the 2SLAQ LRG and QSO sample to make comparisons 
for small and large scale effects in the redshift distortions using 
both the new Scoccimarro expression as well as the Kaiser Limit. 

The small-scale random motions of the galaxies, w ra nd, leads 
to an extension in the 7r-direction of £(<r, n). We denote the magni- 
tude of this extension by (u> 2 ) 1//2 (= o"i2); this is usually expressed 
in a Gaussian form (e.g. dA05) 



/K) = 



2n{w1) 1 / 2 



exp 



2(w|)i/s 



(25) 



Now we can combine these small-scale non-linear z-space distor- 
tions with the Kaiser formulae, and hence the full model for £(<r, n) 
is given by 



£(«r,jr)= / £'[a,TT - w z (l + z)/H(z)]f(w z )dw z 



(26) 



where £'[<7,7T — w z (l + z)/H(z)] is given by equation 1241 and 
f(w z ) by equation [25] Using these expressions and our 2SLAQ 
LRG data, we can calculate f3 and (w 2 ) 1/2 for the LRGs. At this 
juncture, it is important to note the scales we consider in our model. 
As can be seen from the data presented in Section 3, a power-law 
fits the data best on scales from 1 to 20 h~ Mpc. Thus, when 
computing the full model for £(<r, tt) (equation |26t. we only use 
data with 1 < a < 20 /i _1 Mpc and 1 < n < 20 h~ 1 Mpc (as 
shown in Fi gures [TT1 and IT2t . 

Returning to Kaiser (1987), the value of f3 can be used to de- 
termine the bias, b, of the objects in question, 



(27) 



provided you know the values of fl m , where Q m (z) is given by 
fi? n (l + z) 3 



(28) 



for a flat universe. The importance of the bias is that it links the 
visible galaxies to the underlying (dark) matter density fluctuations, 



b 8 n 



(29) 



by e 



where the g and the m subscripts stand for galaxies and mass re- 
spectively. However, the precise way in which galaxies trace the un- 
derlyin g matter distribution is still poorly unde r stood. Recent work 
s.g . lBlanton et al.l (2 006), Sch ulz & Whitel j2006t) . ISmifh et all 
d2007l) and IColes & Erdogdul J2007h suggests that bias is poten- 
tially scale-dependent and we note that we do not take this into 
account in the current analysis. Thus, for our purposes, we restrict 
ourselves to the very simple relation, £ 9 = b 2 £m, where b is the 
linear bias term and leave further investigation of the bias for mas- 



sive galaxies at intermediate redshift to a future paper. On the above 
model assumptions we now proceed to estimate the cosmological 
parameters, fl m and f3. 

4.3 Cosmological Parameters from £ (a, tt) models. 

The ratio of observed angular size to radial size varies with 
cosmology. If we have an object which is known to be isotropic, 
i.e. where transverse and radial intrinsic size are the same, fixing 
the ratio of the intrinsic radial and transverse distances yields 
a relation between the measured radial and transverse distances 
depending on cosmological parameters. This comparison is often 
called the "Alcock-Paczynski" test (Alcock & Paczynski 1979; 
Ballinger, Peacock & Heavens 1996). In order to perform this 
test, we assume galaxy clustering is, on average, isotropic and we 
compare data and model cosmologies. Following H02 and dA05, 
for the following sections, we define several terms. 

(i) The Underlying cosmology - this is the true, underlying, 
unknown cosmology of the Universe. 

(ii) The Assumed cosmology - the cosmology used when 
measuring the two-point correlation function and £(<r, n) from the 
2SLAQ LRG survey. Initially in a redshift survey, the only infor- 
mation available is the object's position on the sky and its redshift. 
In order to convert this into a physical separation, you must assume 
some cosmology. As was mentioned earlier, we have considered 
two Assumed cosmologies, the A (fi m , SIa) = (0.3,0.7) and the 
EdS (fi m , fi A ) = (1.0,0.0) cases. 

(iii) The Test Cosmology - the cosmology used to generate 
the model predictions for £(<r, n) which are then translated into the 
assumed cosmology. 

We compare the geometric distortions in both the data and the 
model relative to the same Assumed cosmology. Thus, the key to 
this technique lies in the fact that when the Test cosmology matches 
the Underlying cosmology, the distortions introduced into the clus- 
tering pattern should be the same in model as in the data. The model 
should then provide a good fit to the data, providing the redshift- 
space distortions have been properly accounted for. We can then 
endeavour to find values of 57 m and /3. We assume that for all fur- 
ther discussions, the cosmologies described are spatially flat and 
choose to fit the variable f2^, hence fixing Q A = 1 — O^. 

The relation between the separations and 7r in the Test and 
Assumed cosmologies (referred to by the subscripts t and a respec- 
tively) is the following (Ballinger et al. 1996, H02, dA05): 

f±O a = — Oa (30) 



7Ti 



/||7Ta 



B a 



(3D 



where A and B are defined as follows (for spatially flat cosmolo 
gies): 

c 1 



A 



dz' 



Ho Jo ^/fiQ +fi<J 1 (l + z') 3 ' 



(32) 



(33) 



In the linear regime, the correlation function in the assumed cos- 
mology will be the same as the correlation function in the test cos- 
mology, given that the separations are scaled appropriately, i.e.: 



£t(°"t, 7T t ) = £a(o" a , 7T a ). 



(34) 
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Details on the fitting procedure are given in dA05 (Section 
7.7). Using this AP-distortion test, we calculate values of O m -/3 
for the assumed A cosmology, and present them in Figure [T3] We 
first note that the constraint here is almost entirely on /3 rather than 
O m . Using the £(r) fit with a (starting) r = 7.45 /i _1 Mpc and 
7 = 1.72, we find that fi m = 0.10±g;f^ and 0(z = 0.55) = 
0.40 ± 0.05 with a velocity dispersion of a — 330kms _1 from a 
\ 2 minimization. We have checked these errors by repeating the 
above calculations on the 32 "jackknife" sub-samples. In order to 
make the jackknife calculations less computationally intensive, the 
velocity dispersion is held fixed at 330 km s _1 in every case. Com- 
paring the error contours in Fig.[T3]with the jackknife estimates, we 
again find that the jackknife errors for /3 at ±0.05 are comparable 
to, if not smaller than, those in the error contours in Fig. Q~3] The 
jackknife error in fi m at ±0.14 is comparable to the error con- 
tour in Fig. [T3] As in Fig. QJJ] this agreement may be surprising 
given that we have ignored the covariance in £(cr, n) points which 
is almost certainly non-negligible. Again we argue that a relatively 
poor % 2 fit at minimum may be responsible, leading to a somewhat 
fortuitous agreement of the formal and jackknife error. But on the 
grounds of the jackknife results we believe that the error contours 
shown in Fig. [T3] are reasonably realistic and we shall quote these 
hear after. 

We have also fitted £(<r, n) assuming an EdS cosmology. In 
principle this should give the same result as assuming the A model. 
We show these fi m — /3 fits in Figure[l4] We find that the best fit is 
now Q m = OAOt°o b 25 and p = Q.45±g;fg (x 2 minimization) with 
a velocity dispersion of a = 330kms _1 . A model with 7 = 1.67 
and a (starting) correlation length of ro = 5.65 /i _1 Mpc is used. 
Thus the /3 and the velocity dispersion values are reasonably con- 
sistent with the previous result. However, the value of fi m assum- 
ing an EdS cosmology, is somewhat higher than the best-fit found 
assuming a A cosmology. We assume that the high degeneracy of 
f2 m coupled with slightly different £(r) models in the two cases is 
causing this discrepancy. The contours in Fig.[l4]certainly suggest 
that the constraint on fi m is much less strict in the EdS assumed 
case. 

We have investigated other systematics in the Q m — j3 fits. 
Returning to an assumed A cosmology, there is some small depen- 
dence on the model assumed for f(r). For example, if the slope 
7 = 1.69 from fitting £(r) in the more limited range 0.4 < 
r < 20h~ 1 Mpc is assumed then we find that fi m = 0.10 ± 0.29 
and f3(z = 0.55) = 0.35 ± 0.16 with a velocity dispersion of 
a — 300kms _1 . Further, if instead of using f(r), w p (a) is used 
with slope 7 = 1.83 over the usual 0.4 < r < 50 /i _1 Mpc 
range, we find that the best-fit model prefers a very low value of 
On = 0.02 ± 0.15 and 0(z = 0.55) = 0.40 ± 0.05 with a veloc- 
ity dispersion of a = 360kms~ 1 . The consistency of these different 
models to give values of f2 m , j3 and a pairwise velocity dispersion, 
albeit at a cost of a very loose constraint on il m , is re-assuring and 
summarised in Table|6] Since w(8) also seems to indicate a flatter 
(7 = -1.67 ± 0.07) slope in the 1 < r < 20 /i _1 Mpc range of 
interest for £(cr, n) we take our 'best bet' estimates to be the values 
for 7 = —1.72 given above. These values also give a good over- 
all fit to £(s). We next introduce a further constraint to break the 
Cl m — /3 degeneracy. 

4.4 Further Constraints on Q'q and f3(z) from LRG 
Clustering Evolution 

iMatsubara & Sutol dl996h and|c room & Shanks! ( fl99^ pointed out 
that by combining low redshift and high redshift clustering infor- 
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Figure 13. Likelihood contours of Q®n-/3(z = 0.55) using the geomet- 
ric method of Alcock-Paczynski test and modelling the redshift-space dis- 
tortions. The best-fit values are Q m = 0.10±° 'fo and P{ z = °-55) = 
0.40 ± 0.05 with a velocity dispersion of a = 330kms —1 . Note how a 
value of Qm ~ 0.3 is not ruled out but also the large degeneracy along 
the Qm direction. A A cosmology is assumed, along with a model where 
7 = 1.72 and a (starting) value of tq = 7.45 h _1 Mpc. 
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Figure 14. Likelihood contours of D. l ^ n -f3(z = 0.55) using the geo- 
metric method of Alcock-Paczynski test and modelling the redshift-space 
distortions, assuming an EdS cosmology. The best-fit values arc — 
°- 4 °ta25 and /3 = 0.45+^u using a model with 7 = 1.67 and a (start- 
ing) correlation length of tq = 5.65 h _1 Mpc. A value of Q m ~ 1.0 lies 
within our 1 a contour but again there is a large degeneracy along the f2 m 
direction. 
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Table 6. Best fitting model values of f2 m , /3 and pairwise velocity dispersion, {w 2 } 1 / 2 , using redshift-space distortions alone and assuming a A cosmology. 



ro 


7 


range / h ^^Mpc 


Measure 




/9 


<U>2)l/2/ kms -l 


7.45 


1.72 


0.4-50 


£(r) 


0.10 


0.40 


330 


7.30 


1.83 


0.4-50 


w p {a) 


0.02 


0.40 


360 


7.60 


1.68 


0.4-20 


CM 


0.10 


0.35 


300 


7.34 


1.80 


0.4-20 




0.10 


0.45 


360 



mation, further constraints on Q m and Q,a would be possible. The 
basic idea described in this section is that the Q Tn :/3(z) degenerate 
set obtained from LRG clustering evolution is different from the 
Q m :/3(z) degenerate set obtained from analysing LRG redshift- 
space distortions; by using these two constraints in combination, 
the degeneracies may be lifted. Thus the way we proceed to break 
the degeneracy is to combine our current 2SLAQ LRG results with 
constraints derived from consideration of LRG clustering evolu- 
tion. 

From the value of the mass correlation function at z = 0, 
linear perturbation theory can be used, assuming a test Q m , to 
compute the value of the mass correlation function in real space 
at z = 0.55. This can then be compared to the measured LRG £(r) 
at z — 0.55 to find the value of the bias b(z = 0.55). The cluster- 
ing of the mass at z — can be determined if the galaxy correla- 
tion function is known, assuming that the bias of the galaxies used, 
b(z — 0), is independent of scale. Fortunately, recent galaxy red- 
shift surveys have obtained precise measurements of the clustering 
of galaxies at z w 0. In practice we shall start from f (s) at z — 
and z = 0.55 and use equation [23] to deri ve £(r) in each case. 

We therefore follow da Angela et al. I d2005h and start by intro- 
ducing the volume averaged two-point correlation function £ where 



J Q S 4irs' 2 ds 



(35) 



We do this so that non-linear effects in the sample should be in- 
significant due to the s 2 weighting, setting the upper limit of the 
integral s = 20 h~ 1 Mpc. To calculate equation [35] at z = 0, we 
use the doubl e-power law form t hat is found by the 2dFGRS to 
describe £ (s) ^Hawkins et alj|2003l . Fig. 6) in the numerator. 

Then, the equivalent averaged correlation function in real- 
space can be determined by 



f(r,* = 0) = 



£(«,« = 0) 



1 + U(z 



0) + i/3(z = 0) 2 



(36) 



where f (s) comes from equation[35]and we take the value of /3 fo r 
the 2dFGRS as (3(z = 0) = 0.49 ± 0.09 dHawkins et ai]|2003h . 
Now the real-space mass correlation is obtained with 



£mass (r, Z = 0) 



£(r,* = 0) 



b(z = 0) 2 ' 

where b(z = 0) is given for each test cosmology by 

Q°J [z = 0) 



b(z = 0) 



(3{z = 0) 



(37) 



(38) 



Once we have determined the real-space correlation function 
of the mass at z = 0, its value at z — 0.55 is obtained using linear 
perturbation theory. Hence, at z — 0.55, the real-space correlation 
function of the mass will be: 



£mass(r, Z = 0.55) 



£mass(r,2 = 0) 

G(z = 0.55) 2 '■ 



(39) 



Here, £ mass (r) is the volume-averaged correlation function (with 
1 < r < 20 /i _1 Mpc) and G(z) is the growth factor of per turba- 
tions, given by linear theory dCarroll et aljl9 92: Peebl esFl980h and 
depends on cosmology - in this case the test cosmology. 

Once the value of £ m ass(r, z — 0.55) is obtained for a given 
test cosmology, the process to find f3(z = 0.55) is similar to the 
one used to find £ m a SS (r, 2 = 0), but now the steps are performed 
in reverse: £(s,z = 0.55) can be measured in a similar way as 
£(s, z — 0). The bias factor at z ~ 0.55 is given by: 



b 2 (z = 0.55) 

where f maBS (r 
tained by: 

£(r,z = 0.55) 



j{r,z = 0.55) 

fmass(r, z = 0.55) ' 

is given by equation [39] and z 

£(s,z = 0.55) 

1 + §/3(z = 0.55) + \p{z = 0.55) 



(40) 



0.55) is ob- 



(41) 



The value of f3{z = 0.55) = ttlf'(z = 0.55)/fe(z = 0.55) will 
be obtained, for the given test value of fl m (z = 0) by solving the 
second order polynomial formed by substitutin g £(r, z = 0.55) 
from equation 40 into equation 41 (see iHovle e t al. 2002). The 
confidence levels on the computed values of /3(z = 0.55) are 
calculated by combining appropriately in quadrature the errors on 
f(s, z = 0.55), f(s, z = 0), /3(z = 0.55) and /3(z = 0). 

Combining this clustering evolution constraint with those 
from z-space distortions breaks the degeneracy in the fl m — f3 plane. 
We can now work out the joint-2 parameter best fitting regions. 
This is shown in Figure [15] where the 1, 2 and 3 sigma error bars 
are plotted (dashed lines). The best fitting 2-parameter calculations 
has fi m = 0.25, j3 = 0.45 denoted by the cross in Figure[T5] When 
we consider the 1 -sigma error on each quantity separately we find, 
Q m = 0.25±o;i°, (3 = 0.45 ± 0.05 with a {w 2 ) 1/2 of 330kms _1 . 
A model £(r) is assumed with 7 = 1.72 and ro = 7.45 ft _1 Mpc, 
as is a A cosmology. 

The case of the combined constraint for the EdS assumed cos- 
mology is shown in FigureQJo] The £(r) model with 7 = 1.67 and 
ro = 5.65 ft _1 Mpc is assumed and we find f2 m = 0.35 ± 0.15 
and (5 — 0.45 ± 0.05. Although the 3-sigma contours still reject 
the EdS model, the rejection is less than in the A assumed case. 
Overall we conclude that the combined constraints on f3 are the 
strongest with /3 = 0.45 ± 0.05 consistently produced whatever 
the assumed cosmology or £(r) model. Though the combined con- 
straints on S7 m are less strong and give fi m w 0.3 ± 0.15, they still 
appear consistent with the standard A model. 

As another check, we can use the ratio to deter- 

mine P from equation [23] (see Figure 1 17b . We assume that j3 is 
scale-independent, the z-space distortions are only affected by the 
large-scale infall and are not contaminated by random peculiar mo- 
tions. Fitting over the scales, 5 < s < 50 /i _1 Mpc, we find 
(3 — 0.47 ± 0.14, which is consistent with our determination us- 
ing the distortions. The 1-ct error comes from a standard x 2 anal- 
ysis using the ratios and their errors; these are derived 
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Figure 15. Joint likelihood contours of Q^-/9(z = 0.55) using the geo- 
metric method of Alcock-Paczynski test, modelling the redshift-space dis- 
tortions and including the evolution of clustering constraints, assuming 
the A cosmology. Here we see that the best-fit joint constraint values are 
fin, = 0.25t°^g, P = 0.45 ± 0.05 (marked with the cross) with a 
(tug) 1 / 2 of SSOkms- 1 . 
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Figure 17. The ratio of the redshift-space correlation function to the real- 
space correlation function, measured from the 2SLAQ LRG survey. We as- 
sume a A cosmology for these measurements and fitting over the scales of 5 
- 50 fi _1 Mpc find that P = 0.47±0.14, in very good agreement with our 
redshift-space distortion/evolution of clustering technique measurements. 




0.2 0.4 0.6 0.8 1.0 

Figure 16. Joint likelihood contours for Q l ^ n -P(z = 0.55) using the geo- 
metric method of Alcock-Paczynski test, modelling the redshift-space dis- 
tortions and including the evolution of clustering constraints, assuming an 
EdS cosmology. The joint best-fit has fi m = 0.35±0.15, P = 0.45±0.05 
(marked with the cross) and a (ui 2 ) 1 / 2 of 330kms — 1 . When the joint con- 
straints are considered, a value of fi m = 1.0 can be ruled out at the 3a 
level. 



from adding the jackknife errors on £(s) and £(r) in quadrature. 
We note that this procedure does not take into account the non- 
independence of the correlation function points, suggesting that the 
relatively large error quoted above on j3 may still be a lower limit. 

The low values of Q m « 0.30 and the value of j3 — 0.45 we 
find from the 2SLAQ LRG survey are in line with what is generally 
expected in the current standard cosmological model. Although the 
constraint on (3 is tight, the constraint on Q m is less so and in partic- 
ular the EdS value is not rejected at 3a when clustering distortions 
only are considered. However, when the combined evolution and 
redshift distortions are considered, the EdS value is rejected at the 
3a level. 

Using equations |40l |4"T1 Q m (z = 0) = 0.30 ± 0.15 and 
/3(z = 0.55) = 0.45±0.05, we find that b(z = 0.55) = fi£ 8 (z = 
0.55)//3(z = 0.55) = 1.66 ± 0.35, showing that the 2SLAQ 
LRGs are highly biased objects. This can be compared with the 
value for SDSS LRGs at redshift z = 0.55 which are found to have 
a value of b = 1.81 ± 0.04 (Padmanabhan et al. 2006, Fig. 13). 
The 2SLAQ LRG value is consistent with this SDSS LRG value; 
of course a slightly lower bias may have been expected for 2SLAQ 
LRGs due to the bluer/lower luminosity selection cut. If we assume 
the value found in recent studies of &7 m (z = 0) = .25 l Cole et al 



20051; [Eisenstein et al.l 120051 ; iTegmark et ail 120061 ; iPercival et al 
2006al lbir then our esimate of b becomes b = 1.56 ± 0.33 



Although we leave discussion about the bias estimate and the 
accuracy of the /3-model to a future paper, at the referee's request, 
we compare the non-linear mass correlation function as num eri- 
cally calculated for the standard cosmology dColm et al]|l999h to 
the 2SLAQ LRG f (r), in Figure|9] The errors in £(er, n) are smaller 
at separations 5 to 20 ft -1 Mpc, than at 1 /i~ 1 Mpc, so our esti- 
mates of bias from £(er, n) are weighted towards these larger scales 
where there appears to be approximate consistency with the relative 
amplitudes of £ ma ss and £(r) in Fig. [9] Thus, as mentioned previ- 
ously, our working assumption from here on will be that there is 
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no effect of scale-dependent bias on our Q m — f3 fits and we leave 
further investigation of this issue for future work. 

Finally, taking the value of b(z = 0.55) = 1.66 ± 0.35, we 
can relate b(z — 0) to 6(0.55) using the bias evolution model dFrvl 
fl996h 

b(z) = 1 + [6(0) - l]G(On(0), f2 A (0), z), (42) 

where G(Qm(0), Qa(0), z) is the linear growth rate of t he den- 
sity perturbations jPeebleslll98blll984ICarroll et al.lll992l) . There 
are many other bias models, but here we are making the simple 
assumptions that galaxies formed at early times and their sub- 
sequent clustering is governed purely by their discrete motion 
within the gravitational potential produced by the matter den- 
sity perturbations. This model would be appropriate, for exam- 
ple, in a "high-peaks" biasing scenario where early-type galax- 
ies formed at a single redshift and their co-moving space den- 
sity then remained constant to the present day. There may be ev- 
idence for such a simple evolutionary history in the observed early- 
type stellar mass/luminosity functio ns (e.g. iMetcalfe et alj|200ll ; 
iBrown et al]|2006l: IWake et alj|2006h . From equation 02] and tak- 
ing 6(0.55) = 1.66, implies a value today of 6(0) = 1.52 at 
Z ~ 0.1. This leads to a predicted correlation length today of 
r (z = 0) = 8.5 ± 1.6 /i _1 Mpc (assuming ACDM) which is 
consistent with the 2dFGRS value of r = 8.0 ± 1.0 /i _1 Mpc 
found from averaging the same two matched luminosity bins from 
Table 2 of iNorberg et al.l <2002f) , and previously used in our Fig. 
7. (But note that the 2dFGRS £(s) shown in Fig. 7 might imply a 
somewhat lower value for the 2dFGRS clustering amplitude in this 
bin than r = 8.0 ± 1.0 /i _1 Mpc.) 

Therefore, these correlation function evolution results suggest 
that there seems to be no inconsistency with the idea that the LRGs 
have a constant co-moving space density, as may be suggested by 
th e luminosity funct ion results. But, we note that the LF results 
of Wa ke et alj J2006I) apply to a colour-cut sample, (where 2SLAQ 
LRGs are carefully matched to SDSS LRGs) whereas our cluster- 
ing results are only approximately matched to the 2dFGRS. It will 
be interesting to see if this results holds when the clustering of the 
exactly matched high and low redshift LRGs are compared (see 
Wake et al., 2007, in prep.). 



5 CONCLUSIONS 

We have performed a detailed analysis of the clustering of 2SLAQ 
LRGs in redshift space as described by the two-point correlation 
function. Our main conclusions are as follows. 

(i) The LRG two-point correlation function, £(s), averaged over 
the redshift range 0.4 < z < 0.8, shows a slope which changes 
as a function of scale, being flatter on small scales and steeper on 
large scales, consistent with the expected effects of redshift-space 
distortions. 

(ii) The best fitting single power-law model to the Real-space 
2-Point correlation function of the 2SLAQ LRG Survey has a clus- 
tering length of ro — 7.45 ± 0.35 /i _1 Mpc and a power-law slope 
of 7 = 1.72 ± 0.06 (assuming a A cosmology) showing LRGs to 
be highly clustered objects. 

(iii) Evidence for a change in the slope of the projected correla- 
tion function, which is a prediction of halo occupation distribution 
(HOD) models, is seen in the 2SLAQ LRG survey results, while a 
stronger feature is observed in the angular correlation function of 
the LRGs. A direct explanation for this remains unclear. 



(iv) From redshift distortion models and the geometric Alcock- 
Paczynski test we find fi m = 0.10to ?o and (3{z = 0.55) = 
0.40 ± 0.05 with a velocity dispersion of a = 330kms _1 , as- 
suming a A cosmology. With EdS as the assumed cosmology, 
n m = 0.40±o;25 and P = 0.45to.?o with the best-fitting velocity 
dispersion remaining at a — 330kms _1 . However, in both cases, 
we also find a degeneracy along the f2 m ass,o-/3 plane. 

(v) By considering the evolution of clustering from z ~ 
to zlrg = 0.55 we can break this degeneracy and find that 
Q m = 0.251^5 and (3 = 0.45 ± 0.05 (with a (™ 2 ) 1/2 of 
330kms _1 ) assuming a A cosmology. When the EdS cosmology 
is assumed, we find Q m = 0.35 ± 0.15 and f3 = 0.45 ± 0.05 
(again (w 2 ) 1//2 = SSOkms" 1 ). When the joint constraints are con- 
sidered, a value of f2 m = 1.0 can be ruled out at the 3a level. We 
believe these estimates of (3(z — 0.55) are reasonably robust but 
the values of f2 m are less well constrained, although the above es- 
timate for S7 m = 0.30 ± 0.15 is in agreement with concordance 
values. 

(vi) If we assume a A cosmology with Q m (z = 0) = 0.3 and 
0(z = 0.55) = 0.45 then the value for the 2SLAQ LRG bias at 
z = 0.55 is 6 = 1.66 ± 0.35, in line with other recent measure- 
ments of LRG bias (Padmanabhan et al. 2006). 

(vii) Assuming this b(z = 0.55) = 1.66 value, and adopt- 
ing a simple "high-peaks" bias prescription which assumes LRGs 
have a constant co-moving space density, we predict ro = 8.5 ± 
1.6 h~ 1 Mpc for LRGs at z ~ 0.1. This is not inconsistent with 
the observed result for luminosity matched 2dFGRS 'LRGs' at this 
redshift. 

The clustering and redshift-space distortion resu lts comple- 
ment the other results from the 2SLAQ Survey e.g. Wa ke et al.1 
d2006h . Wake et al. (2007, in prep.) and da Angela et al. (2006, 
in prep). Luminous Red Galaxies may be considered to be "red 
and dead" but they have recently been realised to be very pow- 
erful tools for both constraining galaxy formation and evolution 
theories as well as cosmological probes. Future projects utilising 
LRGs (e.g. to measure the baryon acoustic oscillations or to study 
LRGs at higher redshift/fainter magnitudes) will give us more in- 
sights into today's greatest astrophysical problems, including the 
epoch of massive galaxy formation and the acceleration of the cos- 
mological expansion. 
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